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INTRODUCTI  ON 

A barotropic  filtered  model (SANBAR)  was  developed  by  Sanders,  et. 
al. , 1975,  and  by  others  for  operational  prediction  of  the  tracks  of  tropical 
storms  at  ranges  out  to  72  hours.  Ill  is  model  has  been  used  since  1968  at 
the  National  Hurricane  Center  (NHC),  where  recent  results  (see  Table  1)^ 
indicate  that  it  performs  competitively  with  other  objective  models  v^ich 
are  credited  (Dunn  et.  al. , 1968)  for  the  slow  improvement  in  the  final  sub- 
jective judgement. 

Table  1 

Homogeneous  sample  of  forecast  position  errors  (nm)  over  period  1973- 


1976. 

Model 

12  hr 

Range 

24  hr 

48 'hr 

72  hr 

NHC  67 

55 

117 

287 

433 

NHC  72 

53 

117 

266 

386 

CLIPER 

55 

123 

268 

369 

NHC  73 

53 

113 

248 

388 

SANBAR 

56 

117 

236 

348 

Number  of  cases 

206 

183 

135 

94 

The  character  of  the  model  was  formed  by  the  belief  that  momentum 
advection  is  the  primary  physical  mechanism  for  motion  of  intense  tropical 
vortices.  Loosely  speaking,  the  assertion  is  that  the  storm  is  "steered" 
by  the  larger- scale  current  in  which  it  is  embedded,  as  suggested  by  Riehl 
and  Haggard  and  Sanborn,  1956,  and  by  Jordon  1952,  among  others.  Hence 
it  seemed  that  the  simple  barotropic  physics  should  be  exploited  before  se- 
rious consideration  was  devoted  to  more  complex  physical  effects. 

1.  From  Director's  memo  on  R and  D activities  at  NHC,  dated  July  15,  1977 


8 


The  failure  of  earlier  attempts  at  barotropic  prediction  (Birchfield,  1960, 
Vanderman,  1962,  Kasahara,  1959,  for  example)  to  achieve  operational  ac- 
ceptance was  regarded  as  due  to  the  difficulty  of  establishing  an  adequate  ini- 
tial large-scale  analysis  on  the  basis  of  rawinsonde-derived  pressure  data 
in  lower  latitudes,  where  errors  are  often  as  large  as  natural  variability. 
Hence  the  SANBAR  model  relies  on  an  analysis  of  wind  observations,  aver- 
aged through  the  depth  of  the  troposphere,  and  makes  no  direct  reference  to 
the  pres  sure -height  data.  Difficulties  with  the  separation  of  the  vortex  from 
the  large-scale  flow  in  the  forecasting  process  (Kasahara,  1959),  led  us  to 
utilize  a relatively  small  150-km  mesh  length  and  to  predict  the  storm  as  an 
integral  part  of  the  total  flow  field. 

Two  problems  had  to  be  dealt  with  immediately:  1)  analysis  over  the 
tropical  oceans  where  rawinsonde  data  are  almost  completely  absent,  and 
2)  assessment  of  the  effect  of  the  storm  circulation,  as  distinct  from  the 
large-scale  influence,  on  soundings  made  in  the  vicinity  of  the  storm  (nec- 
essary for  realistic  construction  of  the  total  initial  flow).  The  first  was 
handled  by  the  provision,  at  first  subjectively  and  later  by  objective  auto  - 
mated  means  (Pike,  1975),  of  "bogus"  wind  data  at  a coarse  array  of  points 
covering  large  portions  of  the  SANBAR  forecast  area.  The  second  problem 
was  first  handled  by  subtracting  from  nearby  wind  observations  a vector 
contribution  from  an  idealized  axi-symmetric  vortex  specified  by  the  geo- 
graphical position  of  its  center,  and  by  its  maximum  wind,  eye  diameter  and 
radius  of  influence.  All  of  these  parameters  except  the  last  are  reasonably 
well  known  initially  in  real  time.  The  radius  of  influence  was  subjectively 
determined,  with  results  that  often  seemed  so  unsatisfactory  that  300nm  was 
adopted  as  an  almost  ubiquitous  nominal  value.  When  this  technique  contin- 
ued to  provide  unreasonable -looking  "residual"  large-scale  winds  from  time 
to  time,  it  was  decidedto  ignore  nearby  wind  soimdings  altogether  and  to  sub- 
stitute, at  the  affected  points  of  the  SANBAR  grid,  first  (Pike,  1972),  the 
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vector  sum  of  the  storm  contribution  described  above  and  a constant  large- 
scale  contribution  ef’''al  to  the  recently  observed  direction  and  speed  of  the 
storm  track,  and  later  (Sanders  et.  a^. , 1975  ) a fixed  stream-function  field 
calculated  from  these  winds  and  the  storm  parameters. 

Aside  from  purely  technical  improvements  in  the  SANBAR  calculations 
two  avenues  seem  open  for  improving  performance.  One  stems  from  the 
improvement  in  the  large-scale  oceanic  data  base  over  the  past  decade,  due 
to  increased  numbers  of  better  wind  observations  from  aircraft  and  especially 
to  large  numbers  of  wind  estimates  now  derived  from  cloud  motions  observed 
by  geosynchronous  satellites.  The  other  road  to  improvement,  however  dif- 
ficult it  has  been  in  the  past,  must  lie  in  the  effective  use  of  wind  observa- 
tions in  the  storm- influenced  region.  This  paper  reports  principally  our 
efforts  along  these  two  lines. 

REGRESSION  ESTIMATION  OF  THE  TROPOSPHERIC  MEAN  WIND 


The  current  database  over  the  oceans  comprises  relatively  dense  cov- 
erage in  the  lowest  2 km,  from  surface  observations  by  ship  smd  from  low 
cloud-motion  observations  by  satellite,  and  in  the  layer  from  9 km  to  12  km, 
from  wind  observations  derived  from  aircraft  navigation  systems  and  from 
high  cloud-motions  observed  by  satellite.  We  must  infer  the  tropospheric 
mean  wind  from  information  in  these  two  layers. 

Thus,  following  preliminary  work  by  Pike  (1975)^  , we  derived  some 
definitive  regression  equations  from  an  extensive  sample  of  data  in  the  NHC 
region  of  forecast  responsibility  during  the  period  June  through  October, 

1971  through  1974.  In  these  equations  rawinsonde  wind  observations  at  850mb 
and  250mb  were  used  to  ^proximate  the  tropospheric  mean  wind  calculated  ^ 

from  the  winds  at  the  ten  mandatory  pressure  levels  from  lOOOmb  to  lOOmb 
in  the  same  soundings.  Our  results,  obtained  from  a total  of  11,682  obser- 
vations in  June  through  October  at  Bermuda,  San  Juan,  Hatteras,  Miami.  j: 

Tampa,  Lake  - 
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Charles,  Brownsville  and  Merida,  are  given  in  Table  2. 


Table  2 


Regression  equations  for  estimated  zonal  and  meridional  components  of 
tropospherically  averaged  winds,  in  NHC  region  of  responsibility. 


^lOOO-lOOmb  - +0-4  + 

' 53u35o 

•3'^"250 

reduction  of 
variance 

.92 

root- me  an 
SQuare  error 

3.4 

A 

^lOOO-lOOmb  ■ 

•45V850 

+ •33v25o 

00 

cn 

3.  1 

Note  that  the  reduction  of  variance  is  substantisd  and  that  the  meridional  com- 
ponent is  somewhat  less  well  represented  (presumably  because  of  noisier 
vertical  structure).  Note  further  that  the  root-mean  square  error,  if  not 
reduced  by  the  analysis  and  initialization  processes,  is  sufficien'i  to  produce 
nearly  all  the  present  state-of-the-art  error  in  a 24-hour  forecast  by  simple 
constant  steering,  even  aside  from  additional  error  introduced  by  use  of  data 
other  than  rawinsonde  observations  in  equations  tailored  for  them.  Evidently 
the  value  of  the  equations  would  be  to  reduce  the  large  error  in  the  occasional 
dreadful  forecast  made  in  near-total  ignorance  far  at  sea. 

In  the  anticipation  that  SANBAR  might  be  used  in  other  regions  of  the 
Northern  Hemisphere,  we  obtained  similar  equations  for  the  Eastern  Pacific 
(from  5,  594  observations  at  Vandenberg  AFB,  Hilo,  Johnston  Island  and 
Midway  Island),  and  for  the  Western  Pacific  (from  10,  145  observations  from 
Guam,  Wake  Island,  Truk,  Ponape,  Kwajalein,  Majure,Yap,  and  Koror).  ^ 
Equations  for  these  two  additional  regions  are  given  in  Table  3.  The  results 
indicate  no  substantial  difference  between  the  Atlantic  and  Eastern  Pacific 
areas.  In  the  Western  Pacific,  however,  the  zonal  equation  is  quite  differ- 
ent, and  it  appears  that  the  vertical  structure  of  the  zonal  component  is 
noisier.  Although  the  reduction  of  variance  in  both  components  is  smaller 

1.  The  former  of  these  sets  may  not  represent  the  wind  structure  in  the  zone 
of  tropical  cyclogenesis,  where  next  to  no  rawinsonde  data  exist,  but  they 
should  be  more  reliable  as  the  storm  approaches  these  populated  locations. 
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Table  3 

Regression  equations  for  Eastern  said  Western  Pacific  regions 


Eastern  Pacific 


reduction  of  root-mean- 
variance  square  error 


^1000-100 

+0.2 

+ 

. 52u85o 

+ . 36u250 

.87 

3.9 

A 

viOOO-100 

+0.  1 

+ 

.46v850 

+ . 36v250 

CO 

00 

3.3 

Western  Pacific 

A 

"lOOO-lOO 

-2.2 

+ 

•■^^850 

' ^^“250 

.77 

3.  3 

II 

o 

0 

1 

o 

o 

o 

.™l 

<> 

-0.4 

+ 

•‘*^''850 

;26v250 

.70 

2.7 

in  this  region,  so  is  the  root -me an -square  error,  indicating  that  the  wind  'S 
less  variable  in  the  Pacific  sample,  but  we  cannot  tell  to  what  extent  the  re- 
duction is  temporal  or  spatial. 

Stratification  of  the  above  samples  into  early-season  (June -August)  said 
late-season  (September-October)  portions  showed  only  minor  differences  in 
the  resulting  regression  equations,  reductions  of  variance  said  root-mesm- 
square  errors. 

The  accuracy  of  estimate  of  the  tropospheric -mean  wind  is,  however, 
sensitive  to  the  number  of  levels  at  which  information  is  available.  We  de- 
rived regression  equations  appropriate  for  the  unfortunate  circumstsaices 
when  only  low-level  or  only  high-level  data  were  available,  said  for  the  op- 
timistic hope  that  wind  information  for  a mid-tropospheric  level  (say,  500mb) 
might  somehow  become  available.  These  equations  are  shown  in  Tsdiles  4 
and  5 in  wdiich  the  data  of  Table  1 are  included  for  comparison.  If  only  one 
level  is  available,  the  mesm  wind  can  probably  be  specified  with  little  or  no 
skill  relative  to  local  climatology,  although  one  is  somewhat  better  off  to  have 
data  in  the  upper  than  in  the  lower  troposphere  . If  the  middle  tropospheric 
data  could  be  added  to  observations  at  the  other  two  levels,  substsuitial  im- 
provement would  be  felt  in  specification  of  the  tropospheric  mesoi  wind. 
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Table  4 


Regression  equations  for  one.  two,  and  three  levels  of  information,  re 
gion  of  NHC  forecast  responsibility. 


One  level; 

reduction  of 
variance 

root-mean 
square  error 

^1000-100^®^°^ 

= +4.8 

■'■•'^^'^850 

. 43 

8.8 

“lOOO-lOO^^®®^ 

1 

11 

+-43U250 

.69 

6.  5 

^1000-100^®®®^ 

= -1. 5 

.34 

6.6 

'^lOOO  -100^^®®^ 

= 0.9 

.60 

5.2 

Two  levels; 


^1000-100<®50-25») 

= +0.4  +.53u35„ 

+. 37U25Q 

.92 

3.4 

i 

*i000-100'85»-250) 

= -0.5  +-'45U55(| 

®®“250 

.85 

3.2 

'f' 

•'I 

Three  levels; 


Table  5 

Regression  equations  for  one,  two,  and  three  levels  of  information  , 
Western  Pacific  region. 

. RV  RMSE 

”1000-100^^^**^  ’ ^ ’^^’*250 

^000-100<”®>‘-®-“'"-*%50  . 24  4.3 

Twq  levels: 

"1000- 100^^°®®'  *-2.2  + . 43Uggjj  + . 32u2gQ  .77  8. 3 

^1000-100^^®®®'^®®^  ”®‘‘^  ^ *^®'^850  *'^® 

Three  levels: 

“l0O0-100<“®®-®“®-"®®''-®-'^-®®“850  ^-®-S00*-=“®250  •®®  ‘•® 

n000-100<«®®-®®®-"®®>”-®-^^-®®''850  ^-®*''500  *-®®''250  •®®  >•’ 

Comparing  Tables  4 and  5,  we  see  large  differences  when  data  are  available 
at  only  one  level,  because  of  the  different  vertical  structure  of  wind  fields 
in  the  two  regions;  but  when  data  are  available  tt  three  levels  the  regression 
equations  and  the  root-mean-square  errors,  are  almost  identical. 

The  reductions  of  variance  suggest  more  skill  than  is  actually  present 
in  the  equations,  because  part  of  the  variance  doubtless  resides  in  differences 
in  the  climatological  average  wind  from  station  to  station  within  each  region. 


14 


This  is  particularly  true  of  the  sample  from  the  Atlantic  sector. 

We  made  no  attempt  to  use  zonal  components  as  predictors  for  meridi- 
onal components  of  the  tropospheric  mean  winds,  or  vice  versa.  Such  an 
attempt  mi^t  yield  useful  results  if  trough  and  ridge  tilts,  for  example,  were 
consistently  northeast-southwest  or  northwest- southeast,  but  substantial  im- 
provement over  what  we  have  already  obtained  seems  unlikely. 

When  satellite  cloud-motion  vectors  are  used  in  place  of  rawinsonde  data 
in  die  two- level-equations  (as  would  be  done  in  practice,  for  example)  we 
estimate  a 50%  increase  in  the  root-mean-square  error  of  specification  of 
the  tropospheric-mean  wind.  The  details  of  this  estimation,  and  of  the  entire 
regression  analysis,  are  given  by  Adams  and  Sanders  (1975). 

SOURCES  OF  ERROR  IN  1975  OPERATIONAL  FORECASTS 

We  undertook  to  study  the  causes  of  large  SANBAR  forecast  errors  in 
the  1975  hurricane  season,  with  the  aim  of  applying  our  regression  equations 
in  revised  predictions.  As  a preamble  to  this  effort,  we  made  revised  fore- 
casts based  on  post-season  "best-track"  initial  positions  and  track  velocities 
( Hebert,  1976)  . As  illustrated  in  Table  6,  these  revised  forecasts  pre- 
sented a substantial  improvement  over  the  original  operational  predictions 
at  ranges  out  to  48  hours.  On  the  other  hand,  the  mean  initial  errors  in 
position  and  track  velocity  (based  on  the  premise  that  the  best-track  infor- 
mation represents  absolute  truth),  suggest  that  it  will  be  extremely  difficult 
to  reduce  the  mean  position  error  in  the  24-hour  forecast  below  75nm,  the 
expected  (or  hoped-for)  error  cited  by  Sanders  and  Burpee,  (1968). 

Incidentally,  the  1975  tracks  were  remarkable  in  two  respects:  only  one 
storm  failed  to  recurve  toward  the  northeast,  and  no  storm  executed  a loop 
or  other  exotic  excursion.  There  was  a slight  tendency  for  operational  pos- 
itions to  lie  westward  of  the  best-track  locationa  and  for  operational  track 
velocities  to  be  insufficiently  northeastward,  both  biases  probably  due  to  the 
forecasters'  reluctance  to  anticipate  fully  the  degreeof  recurvature  and  accel- 
eration which  was  actually  occurring.  In  the  event  of  erratic  storm  tracks. 
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Table  6. 


Comparison  of  operational  and  best-track  forecasts 


Mean  position  error, 
operational  foreca8ts(nm) 

OOhr 

“IT 

12hr 

er 

24hr 

121 

86hr 

l8t 

48hr 

72hr 

393 

Mean  position  error, 
best-track  forecasts (nm) 

0 

80 

89 

182 

224 

376 

Percentage  Improvement 
of  best-track  over  oper- 
ational forecasts. 

100 

25 

18 

16 

14 

4 

Mean  magnitude  of 
error  In  operational 
specification  of  ini- 
tial track  velocity  (kts) 

2.8 

- 

- 

- 

- 

- 

Number  of  forecasts 

74 

87 

58 

51 

44 

33 

operational  errors  in  initial  position  and  track  velocity  would  probably  have 
been  larger. 

Twenty  forecasts  were  chosen  for  revision  on  the  basis  of  reanalysis  of 
the  initial  large-scale  flow  pattern.  These  were  about  equally  split  between 
those  that  were  particularly  good  operationally  and  those  that  were  especially 
bad.  (Were  only  egregiously  poor  forecasts  chosen,  it  is  likely  that  any 
reasonable  alteration  of  method  would  show  improvement,  whether  or  not  it 
had  general  merit. ) Sanders  and  Gordon  (1978)  found  the  large  forecast  err- 
ors to  stem  from  a variety  of  causes. 

One  of  the  cases  analyzed  in  detail,  for  Faye  starting  at  OOOOGMT  , 
September  26,  is  illustrated  in  Figs.  1 and  2.  From  a comparison  of  pre- 
dicted and  observed  tracks  in  Fig.  1,  it  is  seen  first  that  the  slow  predicted 
speed  was  responsible  for  the  large  212-nm  operational  error  at  24  hr,  which 
was  improved  in  the  best-track  prediction  only  by  a more  accurate  specifica- 
tion of  the  initial  track  direction.  Second,  neither  forecast  anticipated  the 
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Fig.  2.  Data  bases  for  oceanic  analysis,  September  26, 
OOOOGMT;  a)  near  surface;  b)  200-mb  level.  Positions  of 
bogus  points  are  shown  by  circled  x's.  Radius  of  influence 
is  shown  by  the  dashed  circle.  Dashed  arrows  indicate 
data  not  received  by  operational  deadline  time. 


dramatic  acceleration  aifter  48  hours,  producing  errors  of  961  nm  and  772 
nm  in  the  operational  and  best-track  predictions  at  72  hours. 

The  unusually  dense  initial  observational  coverage  shown  in  Fig.  2,  com- 
prising mainly  wind  estimates  from  sate  Hite -observed  cloudmotion  vectors, 
precludes  lack  of  data  as  an  explanation  of  the  forecast  errors.  It  appears, 
however,  that  the  numerous  observations  within  the  300-nm  influence  dis- 
tance of  Faye  indicate  a large-scale  flow  toward  the  northwest  at  a speed  in 
excess  of  the  specified  initial  speed.  Application  of  the  regression  equations 
in  Table  1 indicate  in  fact,  a large-scale  speed  of  about  15kt.  in  contrast  to 
the  initial  llkt  specified  initially  in  both  the  operational  and  best-track  pre-  j 

dictions.  The  12-hr  observed  displacement  speed  was  in  fact  17kt.  In  the  | 

present  analysis  procedure,  of  course,  these  wind  data  are  discarded  in  fa- 
vor of  the  specified  initial  speed  . Clearly,  useful  data  are  being  lost. 

The  large  error  at  72  hours  arises  from  another  cause.  Figure  3 shows 
the  initial  large-scale  flow  pattern,  with  its  observed  and  predicted  change. 

The  ridge  which  initially  extended  northwestward  of  the  storm  was  predicted 
to  change  little  during  48  hours,  whereas  in  fact  the  trough  in  the  central 
United  States  advanced  northeastward  to  pick  up  the  accelerating  storm.  The 
forecast  error  was  evidently  not  due  to  fixed  boundary  conditions  in  the 
SANBAR  model,  but  rather  to  the  presence  of  important  baroclinic  effects. 

This  view  is  supported  by  the  portions  of  the  NMC  hemispheric  500-mb  prog- 
nostic charts  shown  in  Fig.  4.  Note  that  the  barotropic  forecast  suffers 
from  the  same  defect  as  the  SANBAR  prognosis,  while  the  baroclinic  PE 
forecast  has  the  right  idea,  as  usual,  but  is  a bit  slow  about  it.  Note  further 
from  Fig.  2 that  the  large-scale  structure  in  the  vicinity  of  the  storm  was 
hardly  barotropic.  The  tropospheric  shear  in  this  case  was  substantial  and 
well  organized,  with  a probable  direct  effect  upon  storm  behavior. 

In  other  instances,  large  SANBAR  forecast  errors  were  found  to  be  at- 
tributable to  failure  to  use  satellite -derived  pres  sure -height  data  poleward 
of  30°N,  to  paucity  of  data  of  all  types,  and  to  fixed  values  of  vorticity  and 
stream  function  on  flie  northern  boundaries  of  the  grid  area. 
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Fig.  3.  Large-scale  initial  flow  pattern,  September  28,  OOOOGMT 
(solid  lines)  and  48-hr  stream-function  changes;  a)  observed  and  b) 
predicted  by  SANBAR . Stream-function  lines  are  at  intervals  of 
3x10  m s“2.  Dashed  lines  indicate  stream-function  rises  and 
dotted  lines  falls,  in  units  of  3xl0®m3s"2. 
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Fig.  4.  NMC  prognostic  charts:  a)  36-hr  barotropic 
forecast  valid  at  OOOOGMT,  September  27;  b)  24-hr  PE 
baroclinic  forecast  valid  at  OOOOGMT,  September  27; 
and  c)  72-hr  PE  forecast  valid  at  OOOOGMT,  Sep  - 
tember  29. 


REVISED  ANALYSES  AND  FORECASTS 
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In  the  sample  of  20  forecasts,  we  have  experimented  with  a revised  an- 
alysis method  designed  to  make  better  use  of  the  growing  satellite -derived 
data  base.  The  idea,  aside  from  using  the  more  reliable  regression  equa- 
tions discussed  above,  is  to  use  ttie  current  statistical  analysis  method 
(Sanders  et.  al. , 1975)  only  >)(iiere  rawinsonde  soundings  are  available. 
Elsewhere,  over  most  of  the  oceanic  regions,  statistical  analysis  via  the 
few  bogus  points  is  to  be  abandoned  in  favor  of  a direct  mean-wind  analysis 
derived  from  current  NHC  analyses  for  the  ATOLL  and  200-mb  levels  (Wise 
and  Simpson,  1971).  These  analyses  are  performed  on  the  SANBAR  compu- 
tational grid  south  of  45°n  and  would  appear  to  be  capable  of  bringing  more 
information  to  bear.  The  regions  of  each  analysis  made  are  shown  in  Fig.  5. 
The  boundary  was  determined  subjectively  after  consideration  of  the  unex- 
plained variance  of  the  statistical  analysis  based  on  rawinsonde  data  and  the 
uncertainty  of  estimation  of  mean  winds  from  two-level  satellite  cloud-motion 
vectors. 

Results  to  date  have  been  disappointing,  since  the  revised  analysis  pro- 
ceuure  has  produced  slightly  larger  forecast  errors  than  in  the  operational 
predictions,  with  highly  similar  forecast  tracks.  On  the  whole,  forecast 
track  directions  are  somewhat  improved,  but  a substantial  slow  bias  has  re- 
appeared in  the  forecasts  (a  problem  which  Sanders  et.^l. , 1976  , deemed 
solved).  The  revised  flow  patterns  are  somewhat  weaker  over  the  oceans. 

We  suspect  coding  errors  may  exist  in  the  revised  program  and  intend  to 
pursue  the  matter  further.  The  change  in  analysis  mode  across  the  bound- 
aries in  Fig.  5 led  to  discontinuities  in  the  wind  analysis.  These  however, 
did  not  appear  to  be  a significant  source  of  forecast  error.  This  work  is 
discussed  in  detail  by  Gordon  (1977). 

INTERPRETATION  OF  STORM  INFLUENCED  WINDS 

It  seems  that  rawinsonde  observations  made  within  the  circulation  of  a 
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tropical  storm,  often  with  great  difficulty  and  at  substantial  hazard  to  the 
observers,  should  be  a valuable  source  of  information  concerning  the  track 
of  the  storm.  Yet  both  the  SANBAR  and  the  Hovermale  MFM  (1975),  dis- 
card all  such  observations.  The  reason,  in  the  case  of  SANBAR,  is  that  we 
have  not  been  able  objectively  to  evaluate  the  effect  of  the  storm  circulation 
with  sufficient  accuracy. 

We  have  recently  explored  the  possibility  of  allowing  such  observations 
to  tell  us  the  "best"  values  of  the  storm  parameters  Xj  , X2  , and  Xg  for 
a symmetric  storm  circulation  , Vg  , of  the  form 

X 

V,  . Xi  {sin[  1) 

In  0 5 

where  Xj  £ , Xg  3 in(r  /300)  ' ^ radial  distance  from 

the  center,  Vmax  maximum  wind  speed  near  the  surface,  K is  the 

proportionality  factor  between  this  wind  and  the  maximum  wind  averaged 
through  the  depth  of  the  troposphere,  and  r^  is  the  radial  distance  of 
from  the  center.  Speeds  are  in  knots,  distances  in  nm,  and  the  maximum 
distance  of  storm  influence  is  300nm. 

We  take  as  the  best  set  of  parameters  that  combination  of  values  which 
yields  the  "smoothest"  set  of  recalculated  (large-scale)  winds  defined  by 
Vr  = Yo  ~ ,3^0  * the  observed  wind.  To  determine  smoothness, 

we  calculate  for  each  observation  within  the  influence  region  of  the  storm, 
the  deviation,  V , of  V_  from  a value,  V,^,  derived  from  planar  fits  of  the 

/X/  /v'  ^ A/  P 

u-  and  v-components  of  the  wind.  The  planes  are  determined  from  the  values 
of^Q  or^j.  at  the  nearest  three  stations  defining  a triangle  enclosing  the 
station  in  question.  (If  the  neighboring  station  is  itself  within  the  influence 
region,  is  used).  The  best  set  of  parameters  is  taken  as  that  which  min- 
imizes the  root- me  an- square  value  of  V over  all  stations  for  each  synoptic 
case. 

Finally,  of  course,  the  values  of  should  provide  a good  specification 
of  the  storm-track  velocity  at  the  time  of  the  observations.  We  would  hope 
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that  this  specification  is  as  accurate  as  the  operational  estimate  made  in 
real  time  by  NHC  for  the  official  advisories.  To  this  end,  a value  of  V 
at  the  location  of  the  storm  center  was  determined,  for  the  optimum  set  of 
storm  parameters,  by  stepwise  screening  regression  for  both  zonal  and 
meridional  wind  components,  given  the  station  values  of^^^  . 

Fifty  data  sets,  for  nine  tropical  storms,  were  chosen  for  study  on  the 
basis  of  the  presence  of  two  or  more  simultaneous  rawinsonde  observations 
within  the  influence  region  of  the  storm.  . Understandably,  these  storms  lay 
within  300nm  or  so  of  the  United  States  coast  and  thus  represented  espe  - 
cially  important  forecast  problems  for  NHC.  The  parameters,  Xj  , Xg  , 
and  X3  were  given  discrete  sets  of  values,  with  resulting  selection  fre- 
quencies as  shown  in  Table  7.  We  note  with  surprise  that  in  about  half  the 
instances  the  implied  value  of  is  no  more  than  35kt,  and  that  the  shape 

of  the  radial  profile  is  very  flat,  as  evidenced  by  small  values  of  X3.  The 
current  operational  SANBAR  value  of  X3  = 1.5  is  exceeded  less  than  20% 
of  the  time  in  the  present  sample.  Study  of  individual  cases  shows  that  the 
tropical  storm  is  often  embedded  in  a relatively  weak  cyclonic  circulation 
of  large  scale,  and  that  really  strong  winds  are  rarely  sampled  by  the  rawin- 
sonde system,  leading  to  these  unexpected  results.  Table  8 confirms  an 
association  of  small  values  of  Xj  and  X3. 

For  each  synoptic  case,  we  compared  the  specified  initial  track  velocity 
emerging  from  the  regression  analysis  with  an  estimate  of  the  actual  initial 
velocity  obtained  from  the  best-track  information.  The  mean  magnitude  of 
the  vector  difference  was  4.  2 knots,  slightly  worse  than  the  probable  dis- 
crepancy between  operationally-specified  and  best-track  initial  velocities, 
as  discussed  earlier.  The  frequency  distribution  of  the  differences  in  Table 
9 shows,  however,  that  operational  accuracy  was  probably  exceeded  about 
half  the  time,  and  that  our  regression  procedure  occasionally  yielded  ex- 
tremely large  errors.  Examination  of  cases  showed  that  these  tended  to  be 
instances  in  which  rawinsonde  observations  were  available  in  only  a single 
quadrant  of  the  storm.  There  was  little  bias  in  specified  direction  or  speed. 


Table  7 


Xj(kt) 

N 

r (nm) 

N 

N 

e 

3 

5 

0 

3 

7 

0.2 

13 

10 

8 

6 

1 

0.4 

5 

15 

6 

9 

5 

0.6 

3 

20 

7 

12 

4 

0.8 

6 

25 

4 

15 

7 

1.0 

2 

30 

1 

18 

1 

1.2 

6 

35 

4 

21 

2 

1.4 

5 

40 

2 

24 

2 

1.6 

2 

45 

3 

27 

2 

1.8 

2 

50 

0 

30 

2 

2.0 

1 

65 

0 

33 

2 

2.2 

0 

60 

0 

36 

1 

2.4 

0 

65 

1 

39 

2 

2.6 

1 

70 

2 

42 

1 

2.8 

0 

75 

3 

45 

1 

3.0 

2 

80 

0 

48 

2 

3.2 

1 

85 

2 

51 

1 

3.4 

0 

90 

5 

54 

1 

3.6 

0 

95 

2 

57 

1 

3.8 

0 

100 

Total 

0 

50 

60 

5 

50 

4.0 

1 

50 

none 
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A deficiency  in  westward  motion  might  be  expected  due  to  our  neglect  of  the 
effects  of  the  latitudinal  variation  of  earth  vorticity,  but  the  effect  is  evidently 
small  enough  to  be  masked  by  other  sources  of  error. 

In  a few  cases  in  which  the  storm  center  was  very  close  to  a sounding 
lication  the  regression  result  was  extremely  sensitive  to  the  position  of  the 
center  and  large  errors  were  likely.  Ten  cases  of  this  type  were  recalcu- 
lated with  the  station  less  than  85nm  of  the  center  excluded,  yielding  much 
improved  results.  Figure  6 shows  such  a case,  in  which  the  error  is  13  knots. 
If  the  storm  center  were  located  17nm  to  the  northnortheast.  however,  and 
the  storm  wind  contribution  to  the  Charleston  observation  were  increased 
from  22  to  28  kts,  then  a perfect  specification  would  have  resulted.  Modest 
asymmetry  in  the  storm  circulation  as  well  as  moderate  position  error  could 
produce  the  large  specification  error.  ^ Thus,  we  are  still  unable  to  make 
constructive  use  of  observations  very  close  to  the  storm  center. 

An  especially  interesting  storm  is  Delia  1973,  which  performed  a loop 
along  the  Texas  Gulf  Coast  before  moving  inland.  The  operational  SANBAR 
24-hour  forecasts  were  very  poor  during  this  time,  for  obvious  reasons.  The 
numerous  observations  within  the  influence  region  were  discarded  in  favor 
of  a straight  uniform  large-scale  flow  representing  the  most  recent  storm- 
track  vector.  Our  new  procedure,  as  illustrated  in  Fig.  7,  shows  excellent 
specification  of  the  track  velocity  during  the  loop  (probably  better  than  the 
velocity  estimated  in  real  time),  as  Delia  and  a larger-cyclone,  denoted  by 
the  heavy  block  letter,  circle  about  each  other.  The  looping  motion  of  Delia 
seems  plainly  accountable  as  a barotropic  process.  The  motion  of  the  larger 
cyclone  is  a moot  question. 

Although  the  new  procedure  has  not  been  incorporated  in  the  SANBAR 
analysis  program,  we  estimated  the  errors  that  would  have  ensued  from  using 
the  specified  track  velocity  as  a 12 -hour  extrapolation  forecast.  The  average 
position  error  was  55  nm,  the  same  as  the  12-hour  errors  shown  in  Table  1, 

1.  A recalculation  with  the  Charleston  observations  excluded  yielded  an 
error  of  three  knots. 


f 


Fig.  6.  Initial  mean  winds  for  July  9,  1959,  OOOOGMT.  Central  posi- 
tion and  influence  region  of  hurricane  symbol  Cindy  are  shown,  respec- 
tively by  the  hurricane  symbol  and  (he  dashed  circle.  Positions  of  the 
storm  twelve  hours  earlier  and  later  are  shown  by  hurricane  symbols. 
The  heavy  arrow  represents  a twelve-hour  linear  displacement  at  the 
specified  velocity.  Plotted  winds  represent  Vp,  also  given  by  numer- 
ical notation.  Values  in  parenthesis  denotes  the  observed  wind,  V.  . 
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Initial  mean  winds  for  hurricane  Della,  1200  GMT,  4 September  1973 
Notation  same  as  in  Pig.  6. 


Fig*  d Initial  mean  winds  for  hurricane  Delia,  OOOOGMT,  6 September  1973 
Notation  S2une  as  in  Fig.  6. 
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Fig.  7 e Initial  mean  winds  for  hurricane  Delia,  1200  GMT,  6 September  1973 
Notation  seune  as  in  Fig.  6. 
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but  probably  better  than  current  capability  in  seasons  with  more  erratic 
tracks  than  1973  to  1976.  This  limitation  aside,  primary  potential  advantage 
of  the  new  procedure  is  that  the  large-scale  flow  in  the  storm- influenced  re- 
gion is  not  constrained  to  be  uniform.  jj 

I 

I 

CONCLUDING  SUMMARY 

For  use  of  satellite  cloud-motion  vectors  in  the  SANBAR  prediction 
model,  we  have  derived  some  definitive  regression  formulas,  for  the  zonal 
and  meridional  components  of  the  wind  averaged  over  the  depth  of  the  trop- 
ical troposphere,  given  data  in  the  lower  or  upper  troposphere  alone,  in  both 
these  layers,  and  in  the  middle  troposphere  as  well  as  in  these  layers.  We 
find  that  data  at  only  one  level  yield  a result  little  better  than  use  of  the  cli- 
matological mean,  while  addition  of  middle-tropospheric  data,  difficult  with 
current  satellite  capability,  would  improve  substantially  upon  estimates  based 
on  lower-and  upper-tropospheric  data.  We  find  some  benefit  in  use  of  separ- 
ate formulas  for  the  Atlantic  and  Pacific  areas,  but  little  in  stratification  into 
earlier  and  later  portions  of  the  tropical-storm  season. 

From  a study  of  SANBAR  forecast  errors,  we  find  that  the  present  model 
suffers  from  inadequate  use  of  observations  within  the  storm- influenced  region, 
from  neglect  of  pressure-height  data  outside  the  tropics,  from  fixed  boundary 
conditions  in  the  middle-latitude  portion  of  the  forecast  grid,  and  from  neglect 
of  baroclinic  effects  in  the  large-scale  flow  pattern  surrounding  the  storm. 

We  recommend  a new  procedure  for  use  when  three  or  more  rawinsonde 
observations  lie  within  the  storm- influenced  region.  This  procedure,  in 
which  the  observations  determine  some  parameters  of  the  storm  circulation 
I itself,  appears  capable  of  specifying  the  initial  storm-track  velocity  about 

. as  well  as  present  subjective  practice,  but  removes  the  present  SANBAR 

assumption  of  a uniform  large-scale  flow  within  the  storm- influenced  region. 

It  should  be  especially  useful  when  erratic  tracks  occur  close  to  landfall. 
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